SUPERWIMPS IN SUPERGRAVITY 



JONATHAN L. FENG 
Department of Physics and Astronomy 
University of California, Irvine, CA 92697, USA 

In supergravity theories, a natural possibility is that neutralinos or sleptons freeze 
ZJi " out at their thermal relic density, but then decay to gravitinos after about a year. 

^ - ^ ' The resulting gravitinos are then superWIMPs — superweakly-interacting mas- 

^rr^ , sive particles that naturally inherit the desired relic density from late decays of 

, conventional WIMPs. SuperWIMP dark matter escapes all conventional searches. 

However, the late decays that produce superWIMPs provide new and promising 
^j^' early universe signatures for cold dark matter. 
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1 Introduction 



One of the many virtues of supergravity theories, such as minimal 
supergravity, 1 is that they naturally provide an excellent cold dark matter 
candidate — the neutralino.^ Neutralinos emerge as the lightest supersym- 
metric particle (LSP) in simple models. Assuming i?-parity conservation, they 
' are stable. In addition, their thermal relic density is naturally in the range 

required for dark matter, and there are promising prospects for both direct 
' and indirect detection. 

, Recently we have explored an a lternative solution to the dark matter 

' problem in supergravity theories .12141 Neutralino dark matter is a possibility 

, only if the gravitino mass TOq is above the neutralino mass. In supergravity, 

Q_i' this need not be the case — and the scalar and gaugino masses are all 

of the order of M^oak ~ {F)/Mpi, where M^cak 100 GcV - 1 TeV, Mpi ~ 
1.2 X 10^^ GeV, and (F) are the weak, Planck, and supersymmetry-breaking 
r~| , scales, respectively. The specific ordering depends on unknown, presumably 

0(1), coefficients. 

The gravitino is therefore the LSP in roughly "half" of parameter space, 
and it may be cold dark matter. Assuming that the universe inflates and then 
^ , reheats to a temperature below ~ 10^ GeV — 10^° GeVjI^the number of grav- 

itinos is negligible after reheating. Then, because the gravitino couples only 
gravitationally with all interactions suppressed by Mpi, it plays no role in the 
thermodynamics of the early universe. The next-to-lightest supersymmetric 
particle (NLSP) therefore freezes out as usual; if it is weakly-interacting, its 
relic density will be near the desired value. However, much later, after 



--#^-10^B-10«s, (1) 
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the WIMP decays to the LSP, converting much of its energy density to grav- 
itinos. Gravitino LSPs may therefore form a significant relic component of 
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our universe, with a relic abundance naturally near nr,M ^ 0.23P Models 
with weak-scale extra dimensions also provide a similar dark matter particle 
in the form of Kaluza-Klein gravitons,-^-^ with Kaluza-Klein gauge bosons or 
leptons playing the role of the decaying wimpP Because such dark matter 
candidates naturally preserve the WIMP relic abundance, but have interac- 
tions that are weaker than weak, we refer the whole class of such particles as 
"superWIMPs." 

The super WIMP possibility differs markedly from previous proposals for 
gravitino dark matter. In the original gravitino dark matter scenario, grav- 
itinos have thermal equilibrium abundances and form warm dark matter!^ 
The required ~ keV mass for such gravitinos is taken as evidence for a new 
intermediate supersymmetry breaking scale {F). Alternatively, weak-scale 
gravitinos may be produced with the correct abundances during reheating, 
provided that the reheat temperature is tuned appropriately. In contrast to 
these scenarios, the properties of super WIMP dark matter are determined by 
the two known mass scales M^cak and Mp\. Super WIMP dark matter there- 
fore preserves the main quantitative virtue of conventional WIMPs, naturally 
connecting the electroweak scale to the observed relic density. In addition, 
the mechanism of gravitino production through late decays implies that the 
super WIMP scenario is highly predictive, and, as we shall see, testable. 



2 SuperWIMP Properties 

As outlined above, superWIMP dark matter is produced in decays WIMP — >• 
SWIMP -I- S, where S denotes one or more standard model particles. The su- 
perWIMP is essentially invisible, and so the observable consequences rely on 
finding signals of S production in the early universe. In principle, the strength 
of these signals depend on what S is and its initial energy distribution. For 
the parameters of greatest interest here, however, S quickly initiates elec- 
tromagnetic or hadronic cascades. As a result, the observable consequences 
depend only on the WIMP's lifetime r and the ave rage total el ectro magnetic 
or hadronic energy released in WIMP decayP^IILlS 13 14 15 16]I3 

In many simple supergravity models, the lightest standard model super- 
partner is a Bino-like neutralino. For pure Binos, 
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This decay width, and all results that follow, includes the contributions from 
couplings to both the spin ±3/2 and ±1/2 gravitino polarizations. These 
must all be included, as they are comparable in models with high-scale su- 
persymmetry breaking. In the limit Am = towimp — "^SWIMP ^ m-swiMP, 
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r{B — > 7G) oc (Am)^ and the decay lifetime is 



t{B 7G) « 2.3 X 10^ s 



100 GeV 



Am 



(3) 



independent of the overall towimp, mswiup mass scale. 

If a slepton is the lightest standard model superpartner, its decay width 
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This expression is valid for any scalar superpartner decaying to a nearly 
massless standard model partner. In particular, it holds for i = e, fi, or 
f, and arbitrary mixtures of the £l and £r gauge cigenstates. In the limit 
Am = mwiMP — JTiswiMP *C mswiMP, the decay lifetime is 

4 



£G) 3.6 X 10^ s 
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1 TeV 
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The electromagnetic energy release is conveniently written in terms of 



CeM = eEM^WIMP , 



(6) 



where eem is the initial electromagnetic energy released in each WIMP decay, 
and IwiMP = nwiMp/n^^ is the number density of WIMPs before they decay, 
normalized to the number density of background photons n^'^ = 2({3)T^ /n^ . 
We define hadronic energy release similarly as Chad = ^hadywiMP- In the 
superWIMP scenario, WIMP velocities are negligible when they decay. We 
will be concerned mainly with the case where 5 is a single nearly massless 
particle, and so we define 
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the potentially visible energy in such cases. 

For the neutralino WIMP case, S = 7. (Possible contributions to hadronic 
decay products are discussed in the last section.) Clearly all of the initial 
photon energy is deposited in an electromagnetic shower, and so 



£EM = , £had — . 



(8) 



For the slepton case, the energy release is flavor-dependent. For the case of 
staus, 
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where the range in sem results from the possible variation in electromagnetic 
energy from tt^ and v decay products. The precise value of Eem is in principle 
calculable once the stau's chirality and mass, and the superWIMP mass, are 
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Figure 1. Predicted values of WIMP lifetime t and electromagnetic energy release Cem = 
eem^wimp in the B (left) and f (right) WIMP scenarios for mgwiMP = 1 GeV, 10 GeV, 
. . . , 100 TeV (top to bottom) and Am = mwiMP — mswiMP = 1 TeV, 100 GeV, . . . , 
100 MeV (left to right). For the f WIMP scenario, we assume e^M = 5^t- 



specified. However, as the possible variation in eem is not great relative to 
other effects, we will simply present results below for the representative value 
of Eem = 2^T- 

The lifetimes and energy releases in the Bino and stau WIMP scenarios 
are given in Fig.^for a range of (towimp. Am). For natural weak-scale values 
of these parameters, the lifetimes and energy releases in the neutralino and 
stau scenarios are similar, with lifetimes of about a year, in accord with the 
rough estimate of Eq. and energy releases of 

Cem - 10~' GeV . (10) 

Such values have testable implications, as we discuss in the following two 
sections. 



3 Big Bang Nucleosynthesis 

3.1 Standard BBN and CMB Baryometry 

Big Bang nucleosynthesis predicts primordial light element abundances in 
terms of one free parameter, the baryon-to-photon ratio 77 = UB/n^- At 
present, the observed D, ""^He, ^He, and ^Li abundances may be accommodated 
for baryon-to-photon ratios in the rang^^ 

r;io = r;/10~^° = 2.6 - 6.2 . (11) 
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In light of the difficulty of making precise theoretical predictions and reducing 
(or even estimating) systematic uncertainties in the observations, this consis- 
tency is a well-known triumph of standard Big Bang cosmology. 

At the same time, given recent and expected advances in precision cosmol- 
ogy, the standard BBN picture merits close scrutiny. Recently, BBN baryom- 
etry has been supplemented by CMB data, which alone yields rjio = 6.1±0.4.EI 
Observations of deuterium absorption features in spectra from hig h redshift 
quasars imply a primordial D fraction of D/H = 2.78jiQ'3g ^ 10~^.U^ Combined 

with standard BBN calculations,!^ this yields 7710 = 5.9 ± 0.5. The remark- 
able agreement between CMB and D baryometers has two new implications for 
scenarios with late-decaying particles. First, assuming there is no fine-tuned 
cancellation of unrelated effects, it prohibits significant entropy production 
between the times of BBN and decoupling. Second, the CMB measurement 
supports determinations of 77 from D, already considered by many to be the 
most reliable BBN baryometer. It suggests that if D and another BBN bary- 
ometer disagree, the "problem" lies with the other light element abundance 
— either its systematic uncertainties have been underestimated, or its value 
is modified by new astrophysics or particle physics. Such disagreements may 
therefore provide specific evidence for late-decaying particles in general, and 
superWIMP dark matter in particular. We address this possibility here. 

In standard BBN, the baryon-to-photon ratio ijiq — 6.0 ± 0.5 favored by 
D and CMB observations predictJ^ 

= 0.2478 ±0.0010 (12) 
^He/H = (1.03 ± 0.06) x 10"^ (13) 
^Li/H = A.ltoi X 10"^" (14) 

at 95% CL, where Yp is the ^He mass fraction. At present the greatest discrep- 
ancy lies in ^Li, where all measurements are below the prediction of Eq. H14f) . 
The ^Li fraction may be determined precisely in very low metallicity stars. 
Three independent studics^^ '^'^ ^find 

^Li/H = l.btoi X 10"'" (95% CL) (15) 
^Li/H = 1.72lgJ8 X 10-1° (1(7 + sys) (16) 
^Li/H = 1.23t°oil X 10^1" (stat + sys, 95% CL) , (17) 

where depletion effects have been estimated and included in the last value. 
Within the published uncertainties, the observations are consistent with each 
other but inconsistent with Eq. (|14|l . with central values lower than predicted 
by a factor of 3 — 4. ^Li may be depleted from its primordial value by astro- 
physical effects, for example, by rotational mixing in stars that brings Lithium 
to the core where it may be burned,!^ but it is controversial whether this 
effect is large enough to reconcile observations with the BBN prediction.!^ 
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Figure 2. The grid gives predicted values of WIMP lifetime t and electromagnetic energy 
release C^em = EEMywiMP in the B (left) and f (right) WIMP scenarios for mgwiMP = 
100 GeV, 300 GeV, 500 GeV, 1 TeV, and 3 TeV (top to bottom) and Am = mwiMP - 
mswiMP = 600 GeV, 400 GeV, 200 GeV, and 100 GeV (left to right). For the f WIMP 
scenario, we assume £em = ^^t- The analysis of BBN constraints by Cyburt, Ellis, 

Fields, and 01iv&i2« excludes the shaded regions. The best fit region with (r, ^em) ~ 
(3 X 10^ s, 10~® GeV), where ^Li is reduced to observed levels by late decays of WIMPs to 
super WIMPs, is given by the circle. 



The other Ug ht e lement abundances are in better agreement. For exam ple, 
a global analysis,Eni using the "high" Yp values of Izotov and Thuan,EHI finds 
— 23.2 for 3 degrees of freedom, where is completely dominated by the 
^Li discrepancy. 

3.2 SuperWIMPs and the ''Li Underabundance 

Given the overall success of BBN, the first implication for new physics is 
that it should not drastically alter any of the light element abundances. This 
requirement restricts the amount of energy released at various times in the 
history of the universe. A recent analysis by Cyburt, Ellis, Fields, and Olive of 
electromagnetic cascades finds that the shaded regions of Fig. [3 are excluded 
by such considerations.^! The various regions are disfavored by the following 
conservative criteria: 

D low : D/H < 1.3 x 10"^ (18) 
D high : D/H > 5.3 x 10"'^ (19) 
■^He low : Fp < 0.227 (20) 
^Li low : ■^Li/H < 0.9 x 10^^° . (21) 
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A subset of superWIMP predictions from Fig. ^ is superimposed on this 
plot. The subset is for weak-scale mswiMP and Am, the most natural values, 
given the independent motivations for new physics at the weak scale. The 
BBN constraint eliminates some of the region predicted by the superWIMP 
scenario, but regions with mwiMP, 'TigwiMP ~ M^reak remain viable. 

The ^Li anomaly discussed above may be taken as evidence for new 
physics, however. To improve the agreement of observations and BBN predic- 
tions, it is necessary to destroy ^Li without harming the concordance between 
CMB and other BBN determinations of 77. This may be accomplished for 
(t,Cem) - (3 X 10^ s, 10-^ GeV)PSThis "best fit" point is marked in Fig.El 
The amount of energy release is determined by the requirement that ^Li be 
reduced to observed levels without being completely destroyed - one cannot 
therefore be too far from the "^Li low" region. In addition, one cannot de- 
stroy or create too much of the other elements. ^He, with a binding threshold 
energy of 19.8 MeV, much higher than Lithium's 2.5 MeV, is not significantly 
destroyed. On the other hand, D is loosely bound, with a binding energy 
of 2.2 MeV. The two primary reactions are D destruction through 7D — > np 
and D creation through 7^IIe DD. These are balanced in the channel of 
Fig. m between the "low D" and "high D" regions, and the requirement that 
the electromagnetic energy that destroys ''Li not disturb the D abundance 
specifies the preferred decay time r ~ 3 x 10^ s. 

Without theoretical guidance, this scenario for resolving the ^Li abun- 
dance is rather fine-tuned: possible decay times and energy releases span tens 
of orders of magnitude, and there is no motivation for the specific range of 
parameters required to resolve BBN discrepancies. In the superWIMP sce- 
nario, however, both r and Cem are specified: the decay time is necessarily 
that of a gravitational decay of a weak-scale mass particle, leading to Eq. (Q, 
and the energy release is determined by the requirement that superWIMPs 
be the dark matter, leading to Eq. (|10|l . Remarkably, these values coincide 
with the best fit values for t and Cem- More quantitatively, we note that 
the grids of predictions for the B and f scenarios given in Fig. El cover the 
best fit region. Current discrepancies in BBN light element abundances may 
therefore be naturally explained by superWIMP dark matter. 

This tentative evidence may be reinforced or disfavored in a number of 
ways. Improvements in the BBN observations discussed above may show if the 
^Li abundance is truly below predictions. In addition, measurements of ^Li/H 
and ^Li/^Li may constrain astrophysical depletion of ''Li and may also provide 
additional evidence for late decaying particles in the best fit regionEmm 
Finally, if the best fit region is indeed realized by WIMP SWIMP decays, 
there a re a number of other testable implications for cosmology and particle 
physicsPElwe discuss one of these in the following section. 
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4 The Cosmic Microwave Background 



The injection of electromagnetic energy may also distort the frequency de- 
pendence of the CMB black body radiation. For the decay times of interest, 
with redshifts z ^ 10^ — 10^, the resulting photons interact efficiently through 
76^ 76^, but photon number is conserved, since double Compton scat- 
tering 7e~ — > 776^ and thermal bremsstrahlung eX — > 6X7, where X is 
an ion, are inefficient. The spectrum therefore relaxes to statistical but not 
thermodynamic equilibrium, resulting in a Bose- Einstein distribution function 



/7(S) 
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^E/{kT}+p, „ I 



(22) 



with chemical potential n 0- 

For the low values of baryon density currently favored, the effects 
of double Compton scattering are more significant than those of thermal 
bremsstrahlung. The value of the chemical potential n may therefore be ap- 
proximated for small energy releases by the analytic expressiorPSI 



^ = 8.0 X 10" 
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In Fig. 131 w e show contours of chemical potential /i. The current bound 
is ^ < 9 X lO-s.EHEIwe see that, although there are at present no indica- 
tions of deviations from black body, current limits are already sensitive to the 
superWIMP scenario, and particularly to regions favored by the BEN consid- 
erations described in Sec. |3| In the future, the Diffuse Microwave Emission 
Survey (DIMES) may improve sensitivities to « 2 x 10"^.^^ DIMES will 
therefore probe further into superWIMP parameter space, and will effectively 
probe all of the favored region where the ^Li underabundance is explained by 
decays to superWIMPs. 



5 Summary and Future Directions 

SuperWIMP dark matter presents a qualitatively new dark matter possibility 
realized in some of the most promising frameworks for new physics. In su- 
pergravity, superWIMP dark matter is realized simply by assuming that the 
gravitino is the LSP. When the NLSP is a weakly-interacting superpartner, 
the gravitino superWIMP naturally inherits the desired dark matter relic den- 
sity. The prime WIMP virtue connecting weak scale physics with the observed 
dark matter density is therefore preserved by superWIMP dark matter. 
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'''wimp (sec) 



''"wimp (sec) 



Figure 3. Contours of /i, parameterizing the distortion of the CMB from a Planckian spec- 
trum, in the (t, ^em) plane. Regions predicted by the superWIMP dark matter scenario, 
and BBN excluded and best fit regions are given as in Fig. |2] 



Because superWIMP dark matter interacts only gravitationally, searches 
for its effects in standard dark matter experiments are hopeless. At the same 
time, this superweak interaction implies that WIMPs decaying to it do so after 
BBN. BBN observations and later observations, such as of the CMB, therefore 
bracket the era of WIMP decays, and provide new signals. SuperWIMP and 
conventional WIMP dark matter therefore have disjoint sets of signatures; we 
have explored the new opportunities for dark matter detection presented by 
super WIMPs. We find that precision cosmology excludes some of the natural 
parameter space, and future improvements in BBN baryometry and probes of 
CMB n distortions will extend this sensitivity. 

We have also found that the decay times and energy releases generic in the 
superWIMP scenario may naturally reduce ^Li abundances to the observed 
levels without sacrificing the agreement between D and CMB baryometry. 
The currently observed ^Li underabundance therefore provides evidence for 
the superWIMP hypothesis. This scenario predicts that more precise BBN 
observations will expose a truly physical underabundance of ^Li. In addition, 
probes of CMB /i distortions at the level of ~ 2 x 10~^ will be sensitive 
to the entire preferred region. An absence of such effects will exclude this 
explanation. 

We have considered here the cases where neutralinos and sleptons decay 
to gravitinos and electromagnetic energy. In the case of sleptons, BBN con- 
straints on electromagnetic cascades provide the dominant bound. For neu- 
tralinos, however, the case is less clear. Neutralinos may produce hadronic en- 
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ergy through two-body decays x ~^ ^G*, /iG, and three-body decays x qqG. 
Detailed BBN studies constraining hadronic energy release may exclude such 
two-body decays, thereby limiting possible neutralino WIMP candidates to 
photinos, or even exclude three-body decays, thereby eliminating the neu- 
tralino WIMP scenario altogether. At present, detailed BBN studies of 
hadronic energy rele ase incorporating the latest data are limited to decay 
times T < 10^ s.^l We strongly encourage detailed studies for later times 
r ~ 10^ s, as these may have a great impact on what superWIMP scenarios 
are viable. 

Finally, the gravitino superWIMP scenario has strong implications for 
the supersymmetric spectrum and collider searches for supersymmetry. In 
particular, the possibility that the lightest observable superpartner is charged, 
long thought to be excluded in supergravity by the problems associated with 
charged dark matter, is in fact viable. Searches at hadron colliders for "exotic" 
long-lived heavy charged particles are therefore probes of supergravity, the 
most conventional of supersymmetric theories. 
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